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ABSTRACT 

We develop a new method of modeling microlensing events based on a Monte Carlo simulation that incor- 
porates both a Galactic model and the constraints imposed by the observed characteristics of the event. The 
method provides an unbiased way to analyze the event especially when parameters are poorly constrained 
by the observed lightcurve. We apply this method to search for planetary companions of the lens in OGLE- 
2003-BLG-423, whose maximum magnification A max = 256 ±43 (or A max = 400 ± 115 from the lightcurve 
data alone) is the highest among single-lens events ever recorded. The method permits us, for the first time, 
to place constraints directly in the planet-mass/projected-physical-separation plane rather than in the mass- 
ratio/Einstein-radius plane as was done previously. For example, Jupiter-mass companions of main-sequence 
stars at 2.5 AU are excluded with 80% efficiency. 

Subject headings: Galaxy: bulge — gravitational lensing — planetary systems — stars: low-mass, brown 
dwarfs 



1. INTRODUCTION 

High-magnification microlensing events are exceptionally 
sensitive to the presence of planetary companions to the lens. 
As the projected separation of the source and the (parent- 
star) lens decreases, the size of the images increases, thus 
enhancing the probability that the planet will pass close 
enough to an image to generate a noticeable deviation in the 
lightcurve (iGould & Loebl 19921) . And if the source gets suf- 
ficiently close to the lens, the lightcurve can be perturbed 
by the central caustic as sociated with the parent star itself 
(Griest & Safizadeh 1998). Groups that monitor microlensing 
events to search for planets are well aware of this enhanced 
sensitivity and so devote special attention to these events. 

It is therefore somewhat surprising that the vast majority 
of events monitored by these groups in the past have not 
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been particularly high-magnification, and the vast majority 
of the observations of the few that did reach high magnifica- 
tion were actually performed well away from the peak, when 
the event had far less sensitivity to planets. In particular, 
of the 43 events monitored by the PLANET (Probing Lens- 
ing Anomalies N ETwork) collaboration ( Al brow et al]l200Tl 
Gaudi ~al 120021) over 5 years, most of the sensitivity to plan- 
ets came from just 5 or 6 events, and most of that fro m the 
near-peak regions of th ese events. iTsapras et alJ ll2003l) and 
ISnoderass et alJ (120041) used the relatively spar se OGLE data 
to put limits on planetary systems, although G audi & Hani 
(2004) have argued that planets could not be reliably detected 
from such data alone. 

The main reason for this apparent discrepancy was simply 
a shortage of microlensing alerts. Hence, at any given time, 
there just were no high-magnification events in progress, or 
at least none near their peak. The available telescope time 
then had to be applied to less favorable events. In addition, 
when devising their observational strategy, PLANET consid- 
ered that they would have to characterize the events they were 
monitoring entirely with their own data. Such characteriza- 
tion is absolutely essential to evaluating the sensitivity of each 
event to planets, and it requires a very large number of obser- 
vations on the wings and at baseline when the event has very 
little sensitivity to planets. 

With the commencement of the Optical Gravitational Lens 
Experiment's OGLE-III project (Udalski et al. 2002), the sit- 
uation is radically changed. Using its dedicated 1.3 m tele- 
scope, large-format (35' x 35') camera, generally excellent 
seeing, and ambitious observing strategy, OGLE-III is alert- 
ing microlensing events toward the Galactic bulge at a rate 
of 500 per season. Since microlensing events are uniformly 
distributed in impact parameter uq, and since peak magnifi- 
cation scales A max ~ Mq 1 (for uq -C 1), this implies that there 
are dozens of events with A max > 10 each year, and a handful 
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with A max > 100. Moreover, OGLE-III photometry is pub- 
licly available (literally hours after it is taken), so there is 
generally no need for microlensing followup groups to mon- 
itor the wings or baseline in order to characterize the event. 
That is, followup observing time can be concentrated on the 
highly sensitive peaks of the high-magnification events. In the 
2003 season, such peaks were occurring almost continuously. 
OGLE-III is therefore generating substantial new opportuni- 
ties for microlensing planet search groups such as PLANET 
jAlbrow et alJll998l). the Microlensing Planet Search (MPS, 
iRhie et alJ 1 9991. 12000I). and th e Microlensing Follow-Up Net- 
work ("uFUN. lYoo et alJ2004T) . 

However, the OGLE-III approach also generates substantial 
challenges. In order to monitor a very large area during the 
2002 and 2003 seasons, OGLE-III returned to each field only 
of order 50 times over the roughly 9 month season. During 
the long nights around 21 June, when the bulge transits near 
midnight, the cadence was relatively high, once every two or 
three nights. But at the edges of bulge season, the rate of 
return dropped as low as once per week or less. Hence, some 
extremely high magnification events appeared quite ordinary 
as of their last observation before peak, and then could be 
recognized for what they were only very close to (or past) 
their peak. Indeed, it is quite possible that their true nature 
as high magnification events could not be recognized at all 
from the OGLE-III lightcurve alone. Thus, without additional 
work, the riches generated by OGLE-III could easily pass by 
unnoticed. Beginning in the 2004 season, OGLE adjusted its 
strategy to concentrate on a reduced number of fields that have 
relatively higher expected event rates. Hence, it is expected 
that the above-mentioned problems will be mitigated in future 
seasons. 

Here we develop a new method of modeling microlensing 
events that incorporates both a Galactic model via a Monte 
Carlo simulation and the constraints imposed by the observed 
characteristics of the event. We apply this method to the 
extreme microlensing event (EME) OGLE-2003-BLG-423, 
which at A max ~ 250, proves to have the highest magnification 
ever recorded among single-lens events. As such, the event 
also has the greatest potential sensitivity to planetary compan- 
ions of the lens, with substantial probability of detecting even 
Neptune mass planets, whose event timescale would typically 
be only about 6 hours. This enhanced sensitivity poses spe- 
cial challenges to the analysis because both the form and am- 
plitude of the impact of such small planets on the lightcurve 
will depend on the relative size of the source compared to 
the Einstein ring. If this relative size were known, it would 
be straightforward to calculate its effect. However, since the 
lightcurve is consistent with a point source, our information 
on the source size is limited. 

Similarly, using the lightcurve data alone the impact pa- 
rameter mo is measured only to about 30%. If uq were known 
much more precisely (as it often is for events with relatively 
bright sources), then we would be able to specify with equal 
precision where in relation to the Einstein ring a planet could 
be and still avoid detection. With our less perfect knowledge 
of uq, however, we must be satisfied with a more probabilistic 
statement about these locations. 

Both of these challenges are likely to be generic to the 
analysis of EMEs. Because such events occur with low 
probability, their sources are most likely to be the relatively 
common main-sequence stars that normally lie unnoticed in 
ground-based bulge images, but which can briefly leap to 
prominence in an EME. Since these main-sequence stars are 



faint and hence small, they will most often avoid finite- 
source effects even in EMEs. Their faintness also induces 
large photometric errors in the wings of the lightcurve, the 
region that must be well-measured to accurately determine 
uq. For similar reasons, these two challenges are likely to 
be key issues in future, even more aggressive, microlens- 
ing experiments tha t aim to detect Earth -mass planets ei- 
ther by space-based JBennett & Rhie 2002) or ground-based 
(Gaudi. Han & Gould 2004) observations. 

In our analysis, we will take as o ur start ing point the method 
pionee red by iGaudi & Sacketll |2000) and Albro wet alJ 
( 2000) , which was then applied to a much larger sample by 
Gau di etalJ 12002). However, we improve upon this method 
in several respects. First, we fix the impact parameter uq at 
a series of different values consistent with the event data and 
evaluate the sensitivity to companions at each uq. To find the 
net sensitivity, we must weight each of these outcomes by the 
relative probability that the actual event had that particular 
uq. Second, we determine these relative probabilities not just 
from the fit to the lightcurve data, but by incorporating the 
results of a Monte Carlo simulation of events toward the ac- 
tual line of sight. For each trial uq, we weight the simulated 
events by how well they reproduce both the observed char- 
acteristics of the lightcurve and the probability that the source 
has the luminosity inferred from the lightcurve combined with 
the Monte Carlo event parameters, as determined from the 
Hipparcos luminosity distribution at the observed color of the 
source. This method not only allows us to more accurately es- 
timate the planetary sensitivity, it also permits us to character- 
ize this sensitivity as a function of planet mass and planet-star 
separation, since each simulated event has a definite lens mass 
(drawn from the adopted mass function) and definite lens and 
source distances (and so definite Einstein radius). In contrast, 
the original approach of Albro wet alJ 12000) yielded sensitiv- 
ities in terms of two lightcurve-fit parameters, the planet-star 
mass ratio and the separation in units of the Einstein radius. 

This method would also permit a similarly rigorous statis- 
tical treatment of finite source effects, since each simulated 
event has a definite ratio of source size to Einstein radius. 
However, based on the Monte Carlo, we show that in the case 
of OGLE-2003-BLG-423, finite-source effects are negligible. 

2. DATA 

OGLE-2003-BLG -423 was alert ed by the Early Warn- 
ing System (EWS, lUdalskil l200l on UT 7:38 14 Sept 
2003, almost exactly 24 hours before the peak on HJD' = 
HJD-2450000 = 2897.8070, and less than 5 hours after the 
triggering observation by the OGLE-III observatory in Las 
Campanas, Chile. While the automated alert did not itself call 
any more attention to this event than the other three that were 
alerted simultaneously, the OGLE web site 10 immediately af- 
fixed a "!" to this event, indicating that it was of special in- 
terest. Moreover, from the data available at the web site one 
could see that the event was alread y 3 mag above baseline and 
rising rapidly. See Albrow (2004) for a Bayesian approach to 
determine whether ascending microlensing events are likely 
to achieve high magnification. 

Immediately following the alert, /iFUN decided to focus 
its observations heavily upon this event. Because the event 
was triggered relatively late in the season when the bulge is 
already west of the meridian at twilight, the time per night that 
it could be observed from any one site was restricted: roughly 

10 http://www.astrouw.edu.pl/~ogle/ogle3/ews/ews.html 
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FIG. 1 . — Lightcurve of microlensing event OGLE-2003-BLG-423 near its 
peak on 15 Sep 2003 (HJD 2452897.8070). Data points with 1 a error bars 
are in I (OGLE: red filled circles; ^tFUN Chile: blue filled triangles; /.iFUN 
Israel: magenta open circles) and V (/iFUN Chile: green filled squares). 
All bands are linearly rescaled so that F s and are the same as the OGLE 
observations, which define the magnitude scale. The solid curve shows the 
best-fit point-source/point-lens model for the /-band curve. 



1.5 hours from the Wise Observatory in Israel and roughly 4 
hours from CTIO at La Serena, Chile. The gap in coverage 
between the two observatories was about 6 hours. Because of 
a communication error, the Chile observations have a gap of 3 
hours the first night, but are then generally spaced at roughly 
1 hour intervals on subsequent nights. 

While OGLE-III normally cycles through many fields (sur- 
vey mode), it can also operate in followup mode when there 
is an event of particular interest. OGLE-2003-BLG-423 was 
immediately designated as such an event, but because of com- 
munication problems, it was not observed the first night fol- 
lowing discovery. However, it was observed one to four times 
per night over the next five nights. 

When combined, observations from these three observa- 
tories provide reasonably good coverage of the peak. See 
Figure ^ The /iFUN data are available at the /iFUN web 
site 11 and the OGLE data are available at the above-mentioned 
OGLE EWS web site. 

From a fit to the first two nights of /iFUN observations, 
it was already clear that the effective timescale of the event 
was very short, t e s = u^e = 0.24 day, where ?e is the Ein- 
stein crossing time and mq is the impact parameter in units of 



1 1 http://www.astronomy.ohio-state.edu/~microfun 



the Einstein radius. Combining this with the Einstein cross- 
ing time ?e = 97 days derived from the OGLE data yielded an 
estimate of A max ~ 1 /uq ~ 400, which would be the highest 
magnification single-lens event ever recorded. Recognizing 
the importance of this event, OGLE and /iFUN worked to- 
gether to develop an observation plan that would allow us to 
characterize it as well as possible. Our principal concern was 
that if OGLE returned to its regular cycle of observations and 
/iFUN stopped observing the event altogether (as both would 
normally do several days after the peak), then the OGLE and 
/iFUN observations might barely overlap in time, meaning 
that the two photometry systems could not be rigidly linked 
into a single lightcurve. To resolve this problem, we agreed 
to each observe the event several times for the next few nights 
(weather permitting) and to both regularly continue observing 
it until it got too close to the Sun. 

There are a total of 278 / band images including 150 from 
OGLE, 78 from /iFUN Chile, and 50 from /iFUN Israel. In 
addition there are 7 V band images from /iFUN Chile, all 
taken near peak for the purpose of determining the color of 
the source. Finally, since /iFUN Chile observations are car- 
ried out with an optical/infrared camera, all V and / images 
from this location are automatically accompanied by H band 
images. However, even at peak, the event was too faint in H 
for these observations to be useful. For each data set, the er- 
rors were rescaled to make x 2 P er degree of freedom for the 
best-fit point-source/point-lens (PSPL) model equal to unity. 
We then eliminated the largest outlier and repeated the process 
until there were no 3 a outliers. This resulted in the elimina- 
tion of 1 OGLE point, 1 /iFUN Chile / point and 1 /iFUN V 
point. In the neighborhood of each of these four outliers, there 
are other data points that agree with the PSPL model, showing 
that the outliers are indeed caused by systematic errors rather 
than revealing unmodeled structure in the lightcurve. The fi- 
nal rescaling factors were 1.13 and 0.82 for OGLE and /iFUN 
Chile /, respectively. The other two observatory/filter com- 
binations did not require renormalization. The descriptions 
of the instruments, observing protoc ol, and reduction proce- 
dures are identical to those given in lYoo et alJ J2004I) . The 
photometry is carried out using the DoPHOT-based PLANET 
pipeline. 

3. POINT-LENS MODELS 

The signature of a planetary companion will usually be a 
brief excursion from an otherwise "normal" point-lens mag- 
nificat ion lightcurve. Indeed, as outlined by Gould & Loebl 
( 1992), it is often possible to estimate the planet's properties 
from the gross characteristics of this deviation. The first step 
to searching for p lanets is therefore to fit the lightcurve to a 
point-lens model (Albrowetal. 2000). However, planetary 
deviations can be strongly affected by the finite size of the 
source, even if the rest of the lightcurve is perfectly consis- 
tent with a point source, which c an lead to degeneracies in 
the interpretation of the deviation (Gaudi & Gould 1997), or 
even to a complete failure to detect the deviation. Hence, we 
begin by presenting the best-fit PSPL model, and then inves- 
tigate to what extent finite-source effects can be detected or 
constrained within the context of point-lens models. 

3.1. Point-Source Point-Lens Model 

We fit the data to PSPL models, defined by three lensing 
geometry parameters (to, uq, and fg) as well as a source flux 
F s and a blended-light flux F\, for each observatory-filter com- 
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Table 1. OGLE-2003-BLG-423 Fit Parameters 
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Table 2. OGLE-2003-BLG-423 Fluxes Parameters (Lightcurve Alone) 
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FIG. 2.— Instrumental CMD of a 6' square around OGLE-2003-BLG-423, 
which has been converted to dereddened magnitude and color by translat- 
ing the centroid of the clump giants (magenta circle) to its known position 
[Iq,(V — /)o]clump = (1.00,14.32). Hipparcos main sequence stars placed at 
Rq = 8 kpc are represented as red points. The source (blue triangle with lcr 
errorbar) is significantly fainter than the Hipparcos stars. 



bination ;. That is, 



F i (t) = F Sti A[u(t)]+F bM A(u): 



u 2 + 2 



(1) 



mVm 2 + 4' 

where [u(t)] 2 = u\ + {t -?o) 2 /?e- The best-fit parameters and 
their errors as determined from the lightcurve data alone are 
shown in Table [2 (also see Fig. 0. In Table [2] we present 
flux-parameters from the lightcurve alone that are rescaled to 
be the same as in the OGLE /-band photometry. We find that 
even though the event is quite long, the source is too faint to 
detect microlensing parallax effects. 

There are several notable features of this fit. First, the im- 
pact parameter is extremely small, uq = 0.00250 ± 0.00072 
implying that the maximum magnification is A max = 400 ± 
115. Second, the source is extremely faint, I s = 22.0 ±0.3. 
The OGLE photometry is not rigorously calibrated, but is be- 
lieved to be accurate to a few tenths. Finally, the errors are 
quite large, roughly 30% for each of ?e, uq, and F s . In fact, 
these errors are extremely correlated: appropriate combina- 
tions of these parameters, f e g = mo^e and F max = F s /uq, have 
much smaller errors, 
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NOTE. — f s is rescaled to be the same as in the OGLE-I photometry 



3.2. Color-Magnitude Diagram 

The first step to understanding the impact of these mea- 
surements and their errors is to place the source on an in- 
strumental color-magnitude diagram (CMD). In Figure [2] we 
have translated the instrumental CMD to place the clump at 
[(y-/)o,/ ] = (1 .00, 14.32), which is t he deredde ned color and 
absolute magnitude of the Hipparcos ( ESA 1997) cl ump when 
place d at the Galactocentric distance, Rq = 8 kpc JYoo et alJ 
2004). The source position (as determined from the model fit) 
is shown as a triangle. Since the source is substantially fainter 
than any of the CMD stars, we also plot the Hipparcos lower 
main sequence in the figure (also placed at Ro). Note that the 
source has a dereddened color (V-I)q = 0.73, almost exactly 
the same as the Sun. The instrumental source color (and so 
the source color relative to the clump) is not model depen- 
dent: it can be derived directly from a regression of V flux 
on / flux, without reference to any model. If the source suf- 
fers similar extinction as the clump, then the source is about 
1.4 mag fainter than the Sun would be if placed at the dis- 
tance to the clump, i.e., Rq. The Sun is somewhat evolved 
off the main sequence, but the source is still more than 1 mag 
fainter than zero-age main-sequence stars of the same color 
and metallicity. This offset between the source and the main- 
sequence (placed at Rq) hardly changes even if one assumes 
a substantial difference between the source reddening and the 
mean reddening toward the clump stars because the reddening 
vector is nearly parallel to the main sequence. 

There are basically only three effects that could contribute 
to this offset. First, the source actually could lie well behind 
the bulge. Second, the source could be relatively metal poor 
and so subluminous compared to solar-neighborhood stars. 
Third, the microlenisng model could be in error, either sta- 
tistically or systematical ly. W e make a more detailed investi- 
gation on this offset in § 14.41 

3.3. Finite-Source Effects 

We now explore a set of models that are constrained to hold 
mo and zo at a fixed grid of values. Here, zo = «o/p* an d 
p* = 0*/6e is the angular size of the source 8* in units of the 
angular Einstein radius 6e- We take account of limb darkening 
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FIG. 3. — Likelihood contours relative to the best-fit point-source/point- 
lens model for various point-lens models with finite-source effects, where 
»o is the impact parameter, zo = «o/p*, and P* is the angular size of the 
source in units of the angular Einstein radius. Both top and bottom panels 
are linear in logzo, but differ in scale. As the source-lens separation increases 
(zo ^1), the contours become independent of zo- Note the isolated minimum 
atzo"~0.8. 



by parameterizing the surface brightness S by, 



So 



= i-r 



l--(l-costf) 



(3) 



where z9 is the angle between the normal to the stellar surface 
and the line of sight. Since the source has almost exactly the 
color of the Sun, we assume solar values for T, 

!>= 0.528, r 7 = 0.368. (4) 

Figure shows A\ 2 contours for the various models plot- 
ted as functions of uq and zo- These contours are essentially 
independent of zo, for zo > 1, i.e., for models in which the 
lens does not pass directly over the source. Note that al- 
though models with zo — 1 are strongly excluded, there are 
some models with zo — 0.8 that are permitted at the 1 a level. 

Even though we cannot rule out the model at the isolated 
minimum (uo — 0.004 and zo — 0.8) based on the lightcurve 
data alone, it is extremely unlikely to describe the event if we 
take account of the observed properties of the event combined 
with a Galactic model (see §0}. For the moment, we therefore 
assume zo > 1 and ignore finite-source effects. Nevertheless, 
as we describe below, the process of recognizing an EME in- 
duces a strong selection bias toward large ?e events and so 
toward those with high Be and/or low relative proper motion 
/i. Hence, in §|4]we will investigate the possibility of zo < 1 
more closely after we evaluate the posterior probability. 

4. MODELING THE EVENT 

We first outline a new method to analyze microlensing 
events that incorporates both a Galactic model via a Monte 
Carlo simulation and the constraints imposed by the observed 
characteristics of the event. This method is completely gen- 
eral and can be applied to any microlensing event. We then 
apply the method to OGLE-2003-BLG-423. 

4.1. General Formalism 



In most scientific experiments, ones seeks to determine the 
posterior probability P{a\ A) of a parameter set a = (ai . . .a„) 
given a data set A. By Bayes' theorem, 



P(a\A) = P rel (A\a)P pri (a), 



(5) 



where P re i(A\a) is the probability of the data given the model 
parameters a and P pr i(a) is the prior probability of the param- 
eters. For microlensing events, 



/> re/ (A|a lc ) = exp[-A x V c )/2], 



(6) 



where a lc is the set of parameters describing the lightcurve and 
Ax 2 (a lc ) is the \ 2 difference relative to the best-fit model. 

Since the prior P pr i(a Xc ) is often assumed to be uniform, 
minimization of x 2 (a lc ) is the usual method to find a best 
parameter set a lc . This procedure is appropriate when the 
lightcurve tightly constrains the parameters, but in general it 
is more correct to take account of the priors. However, the pri- 
ors on some of the lightcurve parameters a lc cannot be directly 
specified. While uq and to can be taken as random variables 
drawn from uniform distributions, ?e is a function of several 
independent physical quantities, namely the lens mass, the 
distances to the lens and source, and the transverse velocities 
of the lens and source. We collectively denote these indepen- 
dent physical parameters as a phys . Therefore, Bayes' theorem 
can be rewritten, 

iV c ,« phys |A) = exp[-A x V c )/2]^ n .(a lc ,a phys ), (7) 

where it is understood that some of the lightcurve parameters 
are determined by the physical parameters. At the end of the 
day, one may be more interested in the physical parameters 
(or some subset of them) than the lightcurve parameters, and 
so after obtaining the general probability distribution given by 
equation Q, one may integrate over the remaining "nuisance 
parameters" to get the probability distribution of a specific 
physical parameter. Indeed, we will do exactly this when we 
evaluate planet sensitivities in § 16.21 



4.2. Relative Likelihood 

To apply this general method to OGLE-2003-BLG-423, we 
first simplify Ax 2 (a lc ). In principle, A\ 2 is a function of 
all five parameters, to,uo,tE,F s , and Ft,. In practice, to is 
extremely well determined from the data, while the remain- 
ing four parameters are all highly correlated. That is, since 
f e ff = «o f E an d F max = Fj uq (and so 7 m ; n ) are very well de- 
termined from the lightcurve data (see eq. 0), their product 
^maxfeff = FJe is also well determined. Moreover, since the 
baseline flux is well determined, F s and Fb are almost per- 
fectly anti-correlated. Hence, once ?e is chosen in a particular 
Monte Carlo realization, / s is also fixed to within 0.008 mag 
and all other parameters are rigidly fixed as well. Therefore, 
the relative likelihood is, 



exp[-A X V c )/2] = exp[-A X 2 (f E )/2], 



(8) 



where Ax 2 (?e) is the \ 2 difference relative to the best-fit 
PSPL model. Since all of the lightcurve parameters are deter- 
mined from the physical parameters via the well-constrained 
lightcurve parameters fo, t e g, and F mdx , Bayes' theorem can be 
rewritten in our case, 



P(a phys |A) = exp[-A x 2 (f E (a phys ))/2] J P /)) ,(a p,iys ) 



physN. 



,phys\ 



(9) 
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FIG. 4. — Microlensing event rate toward the Galactic bulge as a function 
of mass. The upper panel shows event rates for main sequence stars and 
brown dwarfs (thin solid line) and for white-dwarf, neutron-star, and black- 
hole remnants ( dotted line). The total event rate is shown as a thick solid 
line (see Gould 2000). The lower panel shows the posterior probability of 
the Monte Carlo events that takes account of both a Galactic model and the 
Hipparcos-based luminosity distribution at the observed source color. 



4.3. Prior Probability 

In order to estimate the prior probability P pr j(a phys ), we 
Monte Carlo the event, considering all combinations of source 
and lens distances, D[ < D s , uniformly sampled along the 
line of sight toward the source (l,b) = (0 .4961, -5.177 5). We 
choose a lens mass randomly from the iGouldl yOOO) bulge 
mass function, and use Gaussian random variables to assign 
each component of the transverse velocities v± of the lens 
and source. Although this mass function is strictly valid only 
for the bulge, it should be approximately valid for the disk 
as well. This is because at b = -5.2, the line of sight gener- 
ally passes more than a scale height below the Galactic plane, 
where the stars are older (and therefore more bulge-like) than 
they are in the immediate solar neighborhood, where the disk 
mass function is best measured. The event rate for this Monte 
Carlo realization is then, 



T cx p HG (D s )D; P HG(D,)Dfe Ef i, 



(10) 



where the density pug as well as the lens and source velocity 
distribution are as given by the lHan & Gouldl ff99fil 
model. The parameters 



H = \Hs-W\ = 



D„ 



v±,i 
D, 



\AGM ( 1 



1 



c 1 \D, D s 



(11) 



and ?e = 9e//J., are all fixed by the chosen distances, trans- 
verse velocities, and mass. While this Galactic model is not 
a perfect representation of the Galaxy, it is substantially more 
accurate than the uniform prior distribution normally assumed 
in most microlensing analyses. 



We next impose another condition on the prior that con- 
strains a phys . Sin ce ?e is fixed by the Monte Carlo, F s is also 
fixed (see § 14. 21 . By comparing this to the position of the 
clump on the instrumental CMD, one can then determine the 
dereddened flux of the source. Since D s is fixed by the Monte 
Carlo, the absolute magnitude of the source for this Monte 
Carlo realization can also be inferred. The prior probability is 
then proportional to A^Hip, the number of Hipparcos stars with 
this inferred absolute magnitude (and within 0.02 mag of the 
measured source V -I color). Hence, if we restrict attention 
to the k-th realization of the Monte Carlo with physical pa- 
rameters a phys ' k , the prior probability is given by, 



P pn (a phys <V(™ H i P )*- 



(12) 



4.4. Posterior Probability 

Combining equations l|9) and d!2t implies that the poste- 
rior probability that a parameter a, lies in the interval a, £ 
[a;,min,a;,max] is proportional to, 

P {at e [a,-, min ,fl ; , max ]) oc J>(a phys <*), (13) 



(14) 
-a;(a phya )), 



where 

P(a^- k ) = (TNm v )k exp[- A X 2 (f E ,,)/2] 

xe(ai(a ph ^ k )-a Lmin )e(a Ln 

a-i is one of the physical parameters a phys (or possibly a func- 
tion of several physical parameters as would be the case for 
a; G fl lc ), both (rA^Hip)i: and fE.t are implicit functions of a physl , 
and 9 is a step function. 

Letting a, = M, we can evaluate the posterior probability 
distribution for the lens mass. Figure |4] shows both the event 
rate T and the posterior distribution of microlensing events 
toward the Galactic bulge as a function of mass. The over- 
all event rate is shown in the upper panel. The thin solid and 
dotted lines represent events from MS stars and brown dwarfs 
(BDs), and from stellar remnants, respectively. Note that the 
number of objects in the mass function steeply decreases as 
the mass increases (M > 0.7M Q ), and in particular that there 
are no MS stars of M > 1M in the Galactic bu lge because 
such stars have alr eady evolved off jHoltzman et alJ fl998t 
IZoccali et alJl200 0). However, since the cross-section of the 
microlensing event is proportional to M 1 / 2 , remnan ts con- 
tribu te of order 20% of the bulge microlensing events JGouldl 
2000). The posterior distribution for the lens mass is shown in 
the lower panel. Note that high masses are strongly favored, 
possibly because of the long timescale ?e. 

Figure [5] shows the distributions of posterior probabilities 
of various other parameters. The thick and thin solid his- 
tograms represent bulge-bulge and disk-bulge events, respec- 
tively. The upper left panel shows the distribution of im- 
pact parameters (histograms) compared to the distribution de- 
rived from the lightcurve data alone (solid curve). Note that 
the best-fit uq from the lightcurve alone is somewhat lower 
than the peak of the posterior distribution (see §[6ji. The im- 
pact parameter and the apparen t ma gnitude are strongly anti- 
correlated as is discussed in § 13.11 and the lower left panel 
shows the distribution of source dereddened apparent mag- 
nitudes (histograms) compared to the distribution based on 
the lightcurve data alone (as represented by the solid curve in 
Fig-Eland by the position and error bar in Hig - . Taking ac- 
count of the prior probabilities T of the Galactic model and of 
the Mi distribution of Hipparcos stars at the observed source 
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FIG. 5 . — Distributions of Monte Carlo microlensing events toward the OGLE-2003-BLG-423 line of sight. The thick and thin solid lines represent bulge-bulge 
events and bulge-disk events, respectively. The impact parameter, absolute magnitude, dereddened apparent magnitude, and source-lens relative proper motion 
are denoted as uq, Mi,Iq, and /i, respectively, while zo = «o/p*. where p„ is the ratio of the source size to the Einstein radius. The Gaussian curves in the Ho and 
Iq panels represent the probability distributions derived from the lightcurve fit alone, i.e., before applying the constraints from the Galactic model. The dotted 
histograms in the middle panels are the distributions of Hipparcos stars at the color of the source, (V — 7)o = 0.73 (left panel), and lenses in the bulge and disk 
obtained from a Galactic model alone (right panel), respectively. 



color (V-/)o = 0.73 drives the source to somewhat brighter 
apparent magnitudes, but it is still consistent with the result 
derived from the lightcurve data alone. 

The dotted histogram in the middle left panel represents 
the absolute magnitude of Hipparcos MS stars with the same 
color as the source. However, the Monte Carlo events favored 
by the lightcurve are dimmer than the average Hipparcos star 
at /?o (see the middle right panel). The lower right panel shows 
the distribution of zo- As we discuss in §|5] the Monte Carlo 
effectively takes account of the selection effects that push to- 
ward low proper motion (and hence lower zo)- The panel 
shows, however, that the probability that zo is small enough to 
generate significant finite-source effects in a point-lens event 
is extremely small. 

The best-fit lightcurve parameters and their errors are 
shown in Tabled Note that parameters are different at the 2a 
level from those with lightcurve alone, and hence the maxi- 
mum magnification of the event is A max = 256 ± 43. 

5. INFLUENCE OF SELECTION EFFECTS 

As mentioned in §|5J the event was alerted only 24 hours be- 
fore peak. The observation just prior to this triggering obser- 
vation was on 2892.6, about 4 days previous. The event was 
not alerted from that observation because up to that point there 



were only two detections on the subtracted images, whereas 
the alert threshold is set at three in order to avoid spurious 
events. Even had the event been alerted, it would have been 
flagged as having an impact parameter uq = 0.0 ± 0.2 and 
therefore would not have been recognized as a high magni- 
fication event. This 4-day cycle time, which was typical for 
OGLE-III observations in mid-September, introduces signifi- 
cant selection effects in the recovery of EMEs. 

The primary effect is to select for long events. For example, 
if we consider an event with the same source star, same im- 
pact parameter, and same magnification on 2892.6, but with 
?e shorter by a factor 2/3, then it would not have been dis- 
covered until after peak. That is, such an event would also 
not have triggered an alert on 2892.6, but at the 2896.6 ob- 
servation, at which point it would have already been 0.5 days 
past peak. By the time followup observations started, it would 
have been a day past peak, and so magnified only about 40 
times. While still impressive, this would not have garnered 
either the attention or the intensive observations triggered by 
OGLE-2003-BLG-423. 

Thus, the fact that the observed timescale is long com- 
pared to that of typical bulge microlensing events is explained 
largely by selection. However, this selection effect is already 
fully accounted for in the Monte Carlo. Consider a Monte 
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FIG. 6.— Detection efficiency of OGLE-2003-BLG-423 in units of planet- 
star mass ratio q and separation d (normalized to the Einstein radius). The 
upper panel shows the detection efficiency contours (25%, 50%, 75%, and 
95%) by minimizing \~ with respect to fo, ?e, and «o 1 Albrow et a lJ200(D . 
and the lower panel shows the 50% efficiency contours for various fixed «o 
(»o = 0.002 to 0.006). As ug increases, the efficiency decreases monotoni- 
cally. For comparison, we present the 50% contours of the former method 
(solid) and the latter method with uq = 0.002, 0.004 (dashed) in the inset. 



Carlo event that has a t E that is much shorter than the best fit 
in Table 1, say t E = 20 days rather than 97 days. This event 
is assigned a source flux that is lower by a factor ~ 4.5 so 
as to reproduce as well as possible the observed lightcurve. 
In fact, the resulting model lightcurve reproduces the peak 
region extremely well: most of the \ 2 difference comes from 
the post-peak wing, which of course did not enter the selection 
process. Thus, there is no additional selection discrimination 
among the Monte Carlo events. 

The event appears to have several "abnormal" character- 
istics relative to typical events as represented in the Monte 
Carlo, and it is of interest to determine which of these are 
brought about by, or enhanced by selection. The most likely 
source distance is about 2.5 kpc behind the Galactic center. 



For bulge-bulge lensing, there is a general selection effect 
driving toward distant sources because these have larger 9e 
and so larger cross sections. See equation (II 0i . However, as 
shown by the dotted lines in the middle right panel of Figure[5] 
this effect alone pushes the peak of the distribution back only 
1.5 kpc (0.7 kpc for bulge-disk lensing) relative to Ro, not 
2.5 kpc. The pressure toward longer events further selects for 
more distant sources because their larger 6 E make the events 
longer. However, since the FWHM of the prior distribution is 
about 3.5 kpc, the adopted distance would not be extremely 
unlikely in any case. 

Another abnormal characteristic is the faintness of the 
source. Up to a point, the event selection procedure would 
appear to pick out brighter sources. As mentioned above, 
a brighter source would exactly compensate in the selection 
process for a shorter event, and these are more common than 
longer events. However, fainter sources are more common 
than brighter ones, and this is a larger effect. Moreover, as 
the source brightness increases, so does its angular size, and 
this eventually cuts off the peak brightness due to finite source 
effects. Based on Figure|5] however, we have concluded that 
the source is probably nowhere near this threshold, so this 
limitation on source size does not enter as a significant factor. 

Finally, the source appears to be dim for its color. If the 
Hipparcos distribution is representative of bulge stars of so- 
lar color, then this feature would actually be selected against: 
more luminous stars would be both more numerous and, if 
lensed by exactly the same lens, more easily recognized be- 
fore peak. However, it may be that the Hipparcos distribution 
is not representative of the bulge. For example, the stars in 
the outer bulge may have significantly lower metallicity than 
those in the solar neighborhood, and therefore be fainter at 
fixed color as is tr ue of subdwarfs in the solar neighborhood 
(e.g.. lC,oiildl200l . 

We conclude that while a number of the features of this 
event appear unusual at first sight, most are explained in 
whole or in part by selection effects. 

6. SEARCH FOR PLANETS 
6.1. Detection Efficiency 

As discus s ed in § □ iGaudi & Sackettl d2000l) and 
lAlbrow et alJ (120001) have already developed a procedure 
for searching for planets in microlensing lightcurves, and 
Gaudi etall fc002l) have applied this to a sample of 43 events. 
For each event, they considered an ensemble of planetary sys- 
tems characterized by a planet-star mass ratio q, a planet-star 
separation (in units of the Einstein radius) d, and an angle a 
of the source trajectory relative to the planet-star axis. We will 
begin by following this procedure, but will introduce several 
important modifications. 

For a given (d,q), we define the detection efficiency e(d,q) 
as the probability that a companion planetary system de- 
scribed by (d,q) would have produced a lightcurve devi- 
ation inconsistent with the observed OGLE-2003-BLG-423 
lightcurve: 

1 f 2 ^ 

€ ^ d 'l)=^J daO[x 2 (d,q,a)-xlsPL- A xl 1 ], (15) 

where x 2 (d,q,a) is the value of \ 2 evaluated for these three 
parameters, Xpspl i ts best-fit value for the PSPL model, and 
9 is a step function. The (d, q) sampling is 0.1 in the log, 
and the angular step size is set to be Aa = y/q/2 in order to 
avoid missing possible planetary perturbations. We choose a 
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FIG. 7. — Convolved detection efficiency contours of OGLE-2003-BLG-423 as a function of the physical mass and the separation of the planetary companion. 
The efficiency contours from MS+BDs and remnants are shown separately in the left and central panels respectively, while the total is shown in the right panel. 
Contours represent e =25%, 50%, 75% and 95% efficiency. 



conservative threshold, Ax 2 hr = 60 dGaudi et all2002l). In this 
incarnation of the procedure, we follow Albrow et al. (2000) 
and adopt for x 2 (d,q,a) the minimum value of x 2 with these 
three parameters held fixed and allowing all other parameters 
to vary. The results are shown in the upper panel of Fig- 
ure [6] The curves represent detection efficiencies of 25%, 
50%, 75%, and 95%. For q = 0.1, companions with separa- 
tion 0.2 < d < 6 are completely excluded by the data because 
they would produce deviations A\ 2 = x 2 (d,q, a)— Xvspl > 60 
that are not observed. However, the effect of planetary com- 
p anions of ma s s ratio q = 10~ 5 would hardly be discernible. 

iGaudi et all (|2002) discussed a potential shortcoming of 
this approach: if (as in the present case) uq is not well con- 
strained by the data, then it is possible that for the procedure 
to say that certain planetary configurations are permitted by 
the data when in fact they are excluded. For example, suppose 
that the measured impact parameter is uq = 0.00250 ±0.00072 
while the actual value is uq = 0.003. For some value of a, the 
caustic induced by a planet could lie right along the mo = 0.003 
trajectory, but the minimization routine might nevertheless 
find a path that lay 5 a (A^ 2 = 25) from this value at mq = 0.006 
and so avoided the planetary causti c but with Ay 2 < 60 (see 
Fig. 6 and the accompanying text in Gaudi et al. 2002). 

To counter this shortcoming, we evaluate the sensitivity at 
each allowed value of mo- Our search of (d, q, mq) parameter 
space reveals no planets. The best fit is at uq = 0.002, d = 1, 
q = 10~ 3 8 , but the A^ 2 is only -2.9, far short of our adopted 
threshold of Ax 2 hl . = -60. We evaluate the efficiency by mod- 
ifying equation 1151 to become 

1 f 2n 

e(d, q; u ) = — daQ [x 2 {d, q, a; u ) - Xp SPL (mo) - Ax 2 hl ] , 
/7r Jo 

(16) 

where x 2 (d,q,a;uo) and Xpspl( m o) are now evaluated at fixed 
uq. Here, uq is defined as the projected separation of the 
source from the center of the caustic induced by the planetary 
companion. This is the appropriate generalization from the 
point-lens case, in which m is the projected separation from 
the (point-like) caustic at the position of the primarily lens. 
The lower panel of Figure [6] shows 50% contours of detec- 



tion efficiency for several values of uq that are consistent with 
the Monte Carlo simulation (see Fig. Ell. For comp arison, we 
present the 50% contours of the Albrow et al. (2000) method 
(solid) and of our new method with uq = 0.002 and uo = 0.004 
(dashed) in the inset. Althoug h the difference is s mall, we 
find that the previous method of Al brow et alJ ( 1200 0) tends to 
overestimate the detection efficiency. 

While it is comforting that the magnitude of this effect is 
small, its sign is somewhat unsettling. Recall that one moti- 
vation for integrating over uq rather than minimizing with re- 
spect to uo was that under the latter procedure, the trajectory 
could "avoid" planetary caustics and underestimate the sensi- 
tivity. In the present case, however, this effect is outweighed 
by the fact that the most probable value of uq is increased by 
taking into account the Monte Carlo compared to the fit to the 
lightcurve alone. See Figure 5. As discussed in § 1, sensitivity 
to planets generally decreases with increasing uq. 

Note that in this particular case, it is only necessary to inte- 
grate over uo (and not all lightcurve parameters as originally 
envisaged by Gaudi et al. 2002) because once mq is specified, 
a ll the other lightcurve parameters are highly constrained (see 

§E3. 

One might be concerned about finite-source effects during a 
planet-caustic crossing. However, we have repeated the calcu- 
lation including finite-source effects for a variety of (d,q,a) 
combinations and for various plausible source sizes as deter- 
mined from Monte Carlo. We find no significant difference in 
planet detection efficiencies. 

6.2. Constraints on Planets 

Microlensing events provide only degenerate information 
on physical properties of the source and lens except in so far 
as other higher order effects such as finite-source effects and 
parallax are detected. However, our new method based on 
Monte Carlo simulations allows us, for the first time, to break 
the degeneracy and place constraints on planetary compan- 
ions in the planet-mass/physical-separation plane, rather than 
scaling these quantities to the stellar mass and Einstein radius 
as was done previously. 

For a given ensemble of Monte Carlo events with posterior 
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probabilities P k (a phys ), the detection efficiency e can be eval- 
uated as a function of the planet mass m and the planet-star 
projected physical separation r± by, 



e(rj_,m): 



J2 N k=1 e (r ± /d,, k 6 E . k ,m/M k ; UQ , k ) P(g^ k ) 



(17) 



where N is the number of Monte Carlo events and 6e(M, di , d s ) 
is the angular Einstein radius. 

Figure Q shows the resulting detection efficiency with the 
curves representing contours for e =25%, 50%, 75%, and 
95%. The left and middle panels show separate detection effi- 
ciencies for the MS+BDs, and the remnant stars, respectively. 
The total efficiency is shown in the right panel. Since rem- 
nant stars are more massive than MS+BDs, at fixed planet- 
star mass ratio, microlensing events by remnant stars probe 
planets of higher absolute mass. Hence, microlensing is less 
efficient as a probe of planets of remnants than of MS+BDs at 
fixed planetary mass. 

Because our Galactic model favors substantially lower 
blending than is implied by the lightcurve alone, the best 
fit magnitude is reduced from A = 400 to A = 256. Never- 
theless, OGLE-2003-BLG-423 is the highest magnification 
single-lens event recorded to date. Despite this honor, the de- 
tection efficiency is not quite as good as two previous high 
magnification events, MACHO-98-B LG-35 (A mxx ~ 100) and 
OGLE-1999-BUL-35 (A max - 125) jGaudi et alJ 120021 see 



also iBond etalJl200l . This is because our observations do 
not cover the peak of the lightcurve nearly as densely as was 
the case in those two events. Peak coverage is key because 
the perturbations by planets mostly occur during a small time 
interval, basically the full width at h alf maximum around the 
peak of the event ( Rattenburv et al. 2002). 
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